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The phase transition of hard-sphere Heisenberg and Neutral Hard spheres mixture fluids has been 
investigated with the density functional theory in mean-field approximation (MF). The matrix of 
second derivatives of the grand canonical potential Q with respect to the total density, concentration, 
and the magnetization fluctuations has been investigated and diagonalized. The zero of the smallest 
eigenvalue Xs signalizes the phase instability and the related eigenvector x s characterizes this phase 
transition. We find a Curie line where the order parameter is pure magnetization and a mixed 
spinodal where the order parameter is a mixture of total density, concentration, and magnetization. 
Although in the fixed total number density or temperature sections the obtained spinodal diagrams 
are quite similar topology, the predominant phase instabilities are considerable different by analyzing 
x s in density-concentration-magnetization fluctuations space. Furthermore the spinodal diagrams 
in the different fixed concentration are topologically different. 
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I. INTRODUCTION 

The study of the properties of the Heisenberg liquid 
emerged as one of the fascinating theoretical problems 
since it was stimulated by MC simulations of Heisenberg 
fluid [l[ , which strongly hint that the ferromagnetic tran- 
sition of the Heisenberg fluid challenge the traditional 
viewpoint to ferromagnetic transition theory. Nijmei- 
jer et al. found the value of critical exponents of the 
Heisenberg and Ising fluid I2jj differ from the expected re- 
sults for the lattice models. They suggested that the the 
effective exponents is related to Fisher renormalization 
with the fixed density. Recently Mryglod et al. sim- 
ulated the larger Heisenberg fluid system and obtained 
the critical exponents from standard finite size scaling 
theory They argued the effective exponents is cor- 
rected by Fisher renormalization when a thermdynam- 
ics system under a constraint [J]. Although the spin liq- 
uids exhibit interesting phase transition by simulations, 
the formal theory of critical phenomena has not solved 
this problem completely. In fact the spin fluids have a 
complex phase behavior as the coupling between the ad- 
ditional spin degrees and spatial coordinates. Besides 
ordinary gas-liquid phase, such a spin fluid model dis- 
plays paramagnetic gas-ferromagnetic Liquid phase tran- 
sitions, critical end point, and tri-critical point. Many 
models have been proposed such as the discrete Ising 
2 , [a, [q, M , p , continuous XY 3 .[lfj| and Heisenberg fluid 
IS III HI 03 OS EflEj Tavares et al. us- 
ing both a mean field (MF) and a more refined modified 
mean field (MMF) density functional theory have found 
in some regime mixed first-order transition, namely a 
condensation-ordering transition fl2| in the pure Heisen- 
berg fluid. Both theoretical works and MC results showed 



that the first-order transition (i.e., an isotropic vapor 
phase and a ferromagnetic liquid phase) are the mixed 
transition of the ordering fluctuations and density fluc- 
tuations and Li et al. investigated it using the method 
proposed by Chen et al. [H|,[l9[ to characterize the mixed 
phase transition|l7[. 

The phase diagrams of binary spin mixture calculated 
by the mean field theory and Monte Carlo simulation 
show many fascinating phase behaviors which come from 
competition among magnetic, condensation and concen- 
tration fluctuations [20| . So far the microscopically moti- 
vated studies try to comprehend the picture of phase be- 
havior of these fluids. However, as we shall demonstrate 
here, knowledge of the complete phase diagram is essen- 
tial but not enough. In order to understand the existence 
of the mixed transition, we have to classify the phase 
transition characters. So our purpose in the present work 
is to characterize the phase transition along ordcring- 
condensation-demixing transition line. Therefore, we in- 
vestigate a mixture of Heisenberg hard spheres and Neu- 
tral hard spheres (HHNH) using the methodfl8l. Oil to 
characterize phase transition. We address these ques- 
tions using density functional theory in the so-called sim- 
ple mean- field theory. As a result, the phase transitions 
of different character take place in this system, as stud- 
ies of the phase diagram and of the associated critical 
phenomena have been shown. 



II. MODEL 

In this paper we consider a binary mixture fluid of 
Heisenberg Hard-Sphere (HHS) and Neutral Hard-Sphere 
(NHS) with equal diameters a. The pair potential for two 
of such particles at positions ri and r 2 is given by 
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u Q h(ri,wi,r 2 ,w 2 ) = A(a6)w s;5 (r,wi,W2) +u hs (r), (I) 
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where r = |ri ~ r2 1 is the distance between particles at 
positions and r 2 and LUi = 4>i) are the orientations 
of particles with < 9i < tt and < 0, < 2ir. The sub- 
scripts a and b denote the components in fluid considered 
{a(b)= HHS, NHS}, and we define 



A(ab) 



1, a = b = HHS, 

0, a = b = NHS, a^b, 



(2) 



which means all the particles have a diameter yielding 
a repulsive hard sphere potential u hs and only the HHS 
particles attract each other via a Heisenberg potential 
u ss . The repulsive hard sphere interaction and the spin 
part are given by 



u hs (r) 



and 



u ss (r,uji, lo 2 ) = 



where 



J(r) 



+oo, r < a, 
0, r > a. 



0, r<a, 
- J(r)si • s 2 , a < r < r c 
0, r > r c . 



-z(r/cr-l) 



(3) 



(4) 



(5) 



here is a unit vector (|s^ 
spin moment with 



r/a 

= 1) in the direction of the 



Si-S2 = cos(wi,o;2) = cos 6*i cos#2+sin#i sin 6*2 cos(0i— <j> 2 ). 

(6) 

In this paper we set the cutoff distance of the Yukawa 
potential r c for r c — oo and the dimensionless parameter 
z for z = 1. To have ferromagnetic phase favoring parallel 
orientations the coupling constant e here is taken to be 
positive. 



III. DENSITY FUNCTIONAL AND 
MEAN-FIELD APPROXIMATION 

The grand potential free energy Q of a nonuniform spin 
liquid mixtures is the minimum of the functional, 



n [p H (r , w) , PN (r) ,V,T\ = F\p H (r ) w) ) p N (r)] (7) 
+ [ dujdrp H (r,u;){V e H t {r,uj)- p, H ) 

N 



+ J dvp N {v){V e N xt {v) - p N ), 

where index H and N means HHS and NHS, respectively. 
p is the chemical potential and the intrinsic Helmholtz 
free energy of the inhomogeneous fluid F[ph{y, uj), pjv(r)] 



is a unique functional of the densities (r, uj) and pn(t). 
The Helmholtz free energy can be written as |23| 



F\ph(t,w), Pn (t)] ^F hs [p H (r,uj),p N (r)} (8) 

P h (ri , uji)u ss (r, uj x , uj 2 )p H (*2 , u 2 ) 

where F hs [pH{v, uj), PAr(r)] is the Helmholtz free en- 
ergy of Hard-Sphere (HS) mixtures system and 
g(?\ , ui\ , v 2 , ui 2 ; A) is the pair distribution function in a 
system which the particles interact via a pairwise poten- 
tial 



U\(r,ui,U2) = h hs (r) + Xu ss (r, uj\,uj 2 ). 



(9) 



The density of Heisenberg particle P h{y, uj) can be split 
into the number density ps(r) and a normalized factor 
Off(r,u;) as 



p H {r,uj) = p H {r)0 H {r,uj) 



where 



/ 



dwO#(r,w) = 1 



(10) 



(11) 



In the local density approximation the Helmholtz free 
energy of HS mixtures system [22| 

F hs = i Jdrdup H (r,u)[hi(A 3 H 4irp B (T,u))-lll2) 
+ \ J drpN ^ [ ln ( A ^ 4 ^(r)) - 1] 

(3 J (l-r?(r)) 2 

where [3 = is the inverse temperature, A Q is the 

thermal de Broglie wavelength of species a, the total den- 
sity p(r) = ,0^ (r) + p N (r) and ?j(r) = (vr/6)p(r)cr 3 is the 
packing fraction. 

In the mean field approximation where the pair 
distribution function takes its large-distance limit 
g(r%, uji, r 2 , uj 2 ; A) = 1 12], the part of the free energy 
related to the spin-spin interactions in Eq.@ becomes 

F mf = 7} J dr 1 dr 2 dui 1 dui 2 p H (r 1 , uji)u ss (r, uj x , uj 2 )p H (r 2 , ^#3) 

After the decomposition P h{y,uj) = p# (r)0# (r, w), 
the minimum condition of the functional 
£l[pH (r, uj), pn{t)] is equivalent to the simultaneous 
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minimization of the grand canonical functional with 
respect to the number densities 



8tt[p H (r,u>), p N (r),T, n] 

5n{p H (r,uj),p N (r),T, p] 
Sp N {r) 

and the orientational configuration 

5fl[p H (r,uj),p N (r),T, p] 
SO h (v,cj) 



0, 



= 0, 



= 0. 



(14) 



(15) 



In the absence of the external field (Vj£ t (r, u>) = 
Kzt( r ) = 0) the system is homogeneous in position, 
but it could be ordered in orientation. So we have 
Ph(*) = Ph, Pn(t) = Pn and O h {t,w) = H {u). In 
the mean field approximation we obtain the bulk expres- 
sion of the grand-canonical free energy density 

— n MF [0 H (ui), p H> p N ,T, p] = fcs(PH,PN) 

-- Ij m Vff|s| 2 + y J duO H {u)hx[A-KO H {u))- p HPH 

here the Helmholtz free energy density of hard-sphere 
system fc%{PH , Pn) is given by Carnahan and Starling 



fcs(PH,PN)= P -J [H PH K%)-1] 

p At] — 3i] 2 



(17) 



and 



jint _ 



du)Og(ui) s(u;), 

poo 

/ drAnr 2 J{r) = Sirea 3 



(18) 
(19) 



From the equilibrium condition (|15|) and the mean field 
grand-canonical free energy (fT6|) we can obtain the equi- 
librium orientational distribution function 



O h (u)) 



a l3p H J mt s-s 



(20) 



/ duie' 3 P Hj " lts ^' 
Here we restrict to the case Oh(w) — Oh (6) and have 

Oflfw) = O H (0) = — e k ^+ xcos8 , ^ 

47T 



where 



k(x) = ln(ir/ sinhcc), 
x = 0p H J in % 



(22) 
(23) 



with the average magnetization £ = J dioOn (0) cos de- 
termined by 



l3 PH J mt Z 



(24) 



Then the grand function free energy density in an equi- 
librium state is expressed by 

ifiMF [Off(6),pH,p N ,T,fl] = fc S (pH,PN)(25) 
+ \j mt PHe + yKx) - PHPH - PNPN- 



IV. THE SPINODAL OF HEISENBERG HARD 
MIXTURE LIQUIDS 



Ferromagnetic 
Fluctuations 



(16) 
PnPn, 



Sp/p 




Condesation 
Fluctuations 



FIG. 1: (Color Online) A schematic plot of the eigenvec- 
tor which characterizes the phase transition. In the space 
(5p/p,5C,S^) the eigenvector can be described by angles 
<f> and 6 which is in the range 8 e (-90°, 90°] and £ 
( — 180°, 180°]. When 8 — the phase transition is a pure 
ferromagnetic phase transition. The phase transition is pure 
condensation when 8 = 90° and cj> = 0. The phase transition 
is pure demixing when 8 = 90° and <j) = 90°. 

At a stable equilibrium state the grand potential has 
its minimum and its variation with respect to the changes 
of number densities and magnetization should be positive 

M, 



sn/v 



d 2 n 

dpi 



i 

2V 

doJ ipN) 



(W 



d 2 n 



(26) 



2^L6p H 5C 



J N 

d 2 n 



6p N 5£ + 2 



dp H d£, 

d 2 n 



: Sp H Sp 



N 



> 0. 



dp N d£ dpHPNt, 
The variation of the grand potential can be rewritten 
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in a matrix form 

1 / M HH M HN Mjj£ \ ( Sp H 

Sn/V =-(Sp B 5 PN J£ ) M NH M NN M Ni 5p N 

1 \ M m M 6N \h : J \ 5£ 

We want to understand the character of the phase tran- 
sition, which combinations of condensation, phase sepa- 
ration, and ferromagnetic order fluctuations are leading 
to the phase transition. For the following analysis it is 
more convenient to rewrite the matrix in terms of the 
total density p = pn + Pn i the concentration C = Ph/p, 
and the ferromagnetic ordering, £. We reexpress the vari- 
ation of the grand potential in new fluctuations space 

(s P ,sc,so 

( M pp M pC M pi \ ( Sp/p 
Sn/V = - ( Sp/p 6C 6£ ) M Cp M CC M Ci 6C 

2 \ M ip M iC M tt ) \ K 

(28) 

where total one particle density fluctuations 5p = 
5pH + 5pN, the demixing fluctuations SC = S^- = 
P 2 [Pn^Ph — Ph^Pn}- And the matrix M is defined by 

/ M pp M pC M pi \ 
M = M Cp M cc M c? (29) 
V M ip M (C M a / 

which has three eigenvalues Ai,A2,A3 with the corre- 
sponding eigenvectors Xi,X2,X3. The positive eigenval- 
ues Ai, A2, and A3 of a stable state guarantee that the 
grand free energy will increase with respect to any varia- 
tion of total number density, magnetization, and concen- 
tration. If the smallest eigenvalue X s vanishes, the system 
can deviate away from the original state without any in- 
crease of the grand free energy and becomes unstable, 
which is a phase transition. The eigenvector x s corre- 
sponding to the zero ei gen value characterizes this phase 
transition precisely [f8l Il9l| and x s is the order parame- 
ter. Therefore, We investigate the direction of eigenvec- 
tor x s in the (Sp, SC, <5£) fluctuations space, which can 
characterize the type of phase transitions. 

In FiglTj we show a schematic plot of the Euclidean 
eigenvector x s in the space (dp/ p,5C,S£,) with unit vec- 
tors (e x ,e y ,e z ) of the Cartesian system. We normalize 
the eigenvector |x s | = 1 of the zero eigenvalue, and calcu- 
late the angles 9 and <j> of x s in the spherical coordinate 

9 = arccos(x s • e z ), (30) 
6 = arccos ( — X ' s Bx | . (31) 

The angle (j> and 9 describe the portion of the total 
number density, the concentration and the magnetiza- 
tion in the phase transition and is defined in the range 
-90° < 9 < 90°, -180° < < 180°. In general, the 
phase transition of the Heisenberg fluid is a combination 



of condensation and ferromagnetic phase transition. If 
9 = 0, we have a pure ferromagnetic phase transition. 
For 9 near zero and (j> ^ 0, we have a ferromagnetic 
dominant phase transition accompanied by a weak con- 
'densation and demixing fluctuations. If 6 = 90° and 
= 0, we have then a pure gas-liquid phase transition. 
For 9 = 90° and <f> near zero, we have a condensation 
dominant phase transition accompanied by a weak phase 
separation and ferromagnetic phase transition. If 9 = 90° 
and 4> = 90°, we have then a pure demixing phase tran- 
sition. For 9 = 90° and <j> near 90°, we have a demix- 
ing dominant phase transition accompanied by a weak 
ordering and condensation phase transition. Now we in- 
vestigate the phase behavior of the binary fluid by deter- 
mining the border of stable region from the zero point of 
the smallest eigenvalue x s . What we will obtain is the 
so-called spinodal. 

V. RESULTS AND DISCUSSION 




Temperature 

FIG. 2: (Color Online) Plot of the eigenvalues Ai, A2,and 
A3, vs T for the Heisenberg Hard sphere mixture fluid at the 
density p* — 0.8 and the concentration C = 0.8. The Curie 
point and the mixed spinodal point here are on the Curier 
line and the mixed spinodal of Fig. [3jf). 

In FigH]the eigenvalues Ai, A2 and A3 are shown as a 
function of the reduced temperature T* = 1/ ' j3e for the 
reduced density p* = pa 3 = 0.8 and the concentration 
C = 0.8. The eigenvalues A2 and A3 always keep to be 
positive, but the the smallest eigenvalue X s = Ai ap- 
proaches zero when decreasing or increasing the temper- 
ature. The instability at higher temperature is actually 
on the Curie line where 9 = 0. And the smallest eigen- 
value at lower temperature meet the mixed spinodal line 
where 9 = 62° and <f> = 80°. So we have a dominant 
phase separation accompanied by a weaker ordering and 
weakest gas-liquid phase transition. With the zero points 
of A s for different concentrations C we can get the spin- 
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FIG. 3: (Color Online) (a,d) The spinodal curve of the Heisen- 
berg Hard Sphere mixture fluid in different density sections 
(p* == 0.3,0.8). The blue dashed line, the Curies Line; the 
green solid line, the mixed spinodal curve; CEP, critical end 
point; IG, Isotropic Gas; FL, Ferromagnetic Liquid. (b,e) The 
magnetization £ along the mixed spinodal curve with respect 
to concentrations. The pink solid line, £. (c,f) The angles 
along the mixed spinodal curve with respect to concentra- 
tions. The black solid line, the angle 9 ; the red dashed line, 
the angle t/> 



odal phase diagram of the Heisenberg Hard sphere liquid, 
which is shown in Fig|3] 



A. Case: In the isotropic phases £ = 

At the isotropic phase of the Heisenberg fluid, there is 
no total magnetization and x = 0. In this case we have 



lim 

x^0 



d 2 k(x) 
dx 2 



(32) 



For all densities p the element My is positive expect 
the element M^, which can become zero and 9 = cor- 
responds to a pure ferromagnetic phase transition. 

From = we can get the Curie line of the Heisen- 
berg Hard mixture fluids. With the reduced density 
p* = pa 3 and the reduced temperature T* = l/(3e we 
can express the Curie line as 



T = — p* 
3 H 



(33) 



which agrees with the result of Li et al [17j because of 
the simple mean field approximation. And the mean- 
field approximation that makes the phase transition be 



pure ferromagnetic. In a more accurate theory the phase 
transition is not pure ferromagnetic and should be ac- 
companied by the weak condensation and demixing phase 
transition. 



B. Case: £ 7^ 0, the total density is fixed 



When £ 7^ the determinant of the coefficient matrix 



M 



Det 



M pp M pC M pi 
M Cp M cc Ma 
M ip M ic A% 



= 0. 



(34) 



which can be solved numerically. In Fig[3](a), we show 
the Curie line (Blue Dash Line) and the mixed spinodal 
(Green Solid Line) of the HHNH mixtures with the fixed 
density p* = 0.3 section. 

the magnetization and the angle as a function of con- 
centrations. The phase diagram Will be discussed in de- 
tail and characterized the type of phase transitions in the 
method [II]. 

In Fig|3] (a) we obtain the spinodal phase diagram of 
the Heisenberg Hard spheres mixture liquid in the con- 
stant density in the temperature concentration plane. 
The mixed spinodal meets the Curie line at the criti- 
cal endpoint with the concentration C = 0.85 and the 
reduced temperature T* = 2.13. Below the temperature 
of the critical endpoint there is a first-order phase tran- 
sition between isotropic vapor and ferromagnetic liquid. 
The magnetization £ along the mixed spinodal is shown in 
FigJ3](&) by using Eq. (j2"4")) . The more Heisenberg particles 
the the liquid have, it show the bigger the magnetization. 
For the mixed spinodal the phase instability is a combina- 
tion of condensation, demixing and ferromagnetic phase 
transition. To characterize the phase instability precisely 
we investigate the angles <fi and 9 of the eigenvector x s 
along it shown in FigJ3] (c). For concentration near 1.0 
the angle <f> is very small and positive and 9 > 45°, where 
the phase instability is predominantly related to the con- 
densation. The positivity of <f> means that an increase of 
the total number density will be accompanied by a small 
demixing and ordering. This result is plausible, because 
the increase of concentration (more Heisenberg particles) 
will enhances the average attraction of the system which 
results in more ordering and demixing. The magnetiza- 
tion £ for the concentration larger than 0.95 is larger than 
0.5. With the decrease of the concentration the angle <fr 
increases but 9 decreases continuously. When the system 
approaches the critical endpoint C = 0.85, the angle 9 
approaches 0. The phase transition at critical endpoint 
is then a pure ferromagnetic phase transition in the mean 
field approximation. 

The spinodal curve and the CEP with C = 0.69 and 
T* = 2.89 at p* = 0.8 section is displayed in Figg] (d). 
Corresponding the magnetization in Fig|3](e) and the an- 
gles in Fig[3](/) have a obvious region 0.735 < C < 0.965 
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that 8 > 45° and <f> > 45°, in which the phase instability 
is predominantly related to the demixing phase transi- 
tion accompanied by the ordering and condensation. By 
comparing Fig[3] (c) with (/) the demixing predominant 
region disappear in the lower total density section. 



sition predominate in the higher temperature section. 
D. Case: £ 7^ and Concentration is fixed 



C. Case: £ 7^ and Temperature is fixed 
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FIG. 4: (Color Online) (a,c) The spinodal curve of the Heisen- 
berg Hard Sphere mixture fluid in different Temperature sec- 
tions (T = 4.5,2.5). The blue dashed line, the Curies Line; 
the green solid line, the mixed spinodal curve; CEP, criti- 
cal end point; IG, Isotropic Gas; FL, Ferromagnetic Liquid. 
(b,d) The angles along the mixed spinodal curve with respect 
to the total density. The black solid line, the angle 9 ; the red 
dashed line, the angle <j>. 

The spinodal diagrams in the density concentration 
plane are shown in Fig. |Ua,c). In order to elucidate 
the differences of the mixture phase behavior at different 
temperatures, phase behavior is characterized by 8 and <p. 
Fig. [3Jh) show in the T — 4.5 section the angles 4> > 85° 
and 8 < 45° along the mixed spinodal curve, in which the 
phase instability is predominantly related to the ordering 
phase transition accompanied by the weak demixing and 
tiny condensation. But in the lower temperature section 
T = 2.5 ( Fig0Jd) ) the instability predominantly related 
to the demixing accompanied by the weak ordering and 
gas-liquid phase transition at p* > 0.6. Although the 
different temperature sections are quite similar topology, 
the predominant phase instability is considerable discrep- 
ancy between them. Recently, FenzeZ at. presented con- 
centration pressure phase diagrams of the ideal Ising mix- 
ture fluids at different temperatures via the Gibbs ensem- 
ble Monte Carlo simulation [2(| , which the phase behav- 
ior can qualitatively agree with the T sections calculated 
in our model. On the other hand the demixing phase 
transition in the mixtures can be easily observed in the 
lower temperature section but ferromagnetic phase tran- 
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FIG. 5: (Color Online) (a,d,g) The spinodal curve of the 
Heisenberg Hard Sphere mixture fluid in different Concentra- 
tion sections (C = 0.8,0.9,0.99). The blue dashed line, the 
Curies Line; the green solid line, the mixed spinodal curve; 
CEP, critical end point; IG, Isotropic Gas; FL, Ferromagnetic 
Liquid. (b,e,h) The magnetization £ along the mixed spinodal 
curve with respect to the total density. The pink solid line, £. 
(c,f,i) The angles along the mixed spinodal curve with respect 
to the total density. The black solid line, the angle 9 ; the red 
dashed line, the angle <j). 

Density-temperature spinodal diagrams are plotted in 
Fig.([5])(a,d,g) where the different sections belong to dif- 
ferent values of the parameter C. In Figj5] (a) the crit- 
ical endpoint is at p* = 0.345 and T* = 2.28 and the 
slope of the mixed spinodal is positive which was found in 
asymmetric binary dipolar mixtures 2l|. In Fig[5](c) the 
phase instability is predominantly related to the demix- 
ing in the region p* > 0.55 and the ordering in the range 
of 0.345 < p* < 0.55. But the condensation fluctua- 
tions is always weakest in the C — 0.8 section. The 
positive slope of the mixed spinodal means that increase 
of the total number density, at fixed concentration, the 
first order phase transition temperature increase with in- 
creasing p* . When considering mixtures at C = 0.9 one 
finds from Fig. ©(d) that the slope of the mixed spin- 
odal move towards flat. Since increase of the Heisenberg 
particle concentration implies that the condensation in- 
stability gradually predominate the phase transition, as 
shown Fig. ©(f). In the limit C — > 1.0 one recovers the 
MF spinodal diagram of the pure Heisenberg fluid con- 
tain only two fluid phases, which is an isotropic vapor 
and a ferromagnetic liauid[T7j. Here we set C = 0.99, 
and then the spinodal diagram is shown in FigJS] (g) and 
the angle <f> is always near the zero and 8 is large than 45° 
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when p* > 0.25 in Fig[5] (i). The predominant instabil- 
ity is related to the condensation with small ordering and 
demixing fluctuations which coincide with results of pure 
Heisenberg fluid. In particular, the angle is obtained as 
a function of the concentration in Fig[3] thereby provid- 
ing an estimate of the crossover between the demixing 
predominant and the condensation predominant phase 
transition as expected. Such the crossover phe nomena 
can be precisely characterized by the method |l9|. 



imation should be corroborated Monte Carlo simulations 
and further explore the combinations of phase transitions 
on the Curie line with more accurate approximation. 
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VI. CONCLUSIONS 

In this paper we have investigated the phase transition 
of HHNH mixtures by using the density functional theory 
in the mean-field approximation. The phase instability 
of the system is discussed with the method developed in 
Ref . ( [191] ) . From the matrix of the second derivatives of 
the grand canonical free energy f2 with respect to the 
total particle density, the concentration and the magne- 
tization, we can determine the thermodynamic stability 
of the system. When the smallest eigenvalue of the ma- 
trix becomes zero, the system becomes unstable. The 
eigenvector corresponding the zero eigenvalue can char- 
acterize the phase instability precisely. In the total den- 
sity, concentration, magnetization space (5p/p,5C,6t;) 
the normalized eigenvector can be described by two an- 
gles 9 and <j> which is in the range —90° < 9 < 90° and 
-180° < (j) < 180°, respectively. 

For temperature above of the critical endpoints, the 
angle 9 = and the phase transition is pure ferromag- 
netic. This is the result in the mean-field approximation. 
With a more accurate approximation the angle 9 does 
not equal exactly to and is just near 0. Then the total 
particle density and concentration will be related to the 
phase transition, which is predominantly ferromagnetic 
and accompanied by a weak condensation and demixing. 
Below CEP there is a first-order phase transition between 
isotropic gas and ferromagnetic liquid. The phase insta- 
bility along the spinodal near the ferromagnetic liquid is 
a combination of condensation, demixing and ferromag- 
netic phase transition. The phase behavior can topologi- 
cally be worthy to be compared with the simulations [20j 
in the different T sections. Although the different tem- 
perature sections are similar topology, in the lower tem- 
perature section T = 2.5 the instability predominantly is 
related to the demixing but in the T = 4.5 section the or- 
dering phase transition is predominant. The slope of the 
mixed spinodal is the positive in the less concentration 
C = 0.8 sections and is negative in the C = 0.99 section. 
Hence by accurately investigating the angles we quanti- 
tatively describe that the phase instability continuously 
changes from predominant demixing phase transition to 
predominant gas-liquid phase transition. Although these 
crossover phenomena in the concentration sections may 
be considered as different phase transitions, our analysis 
shows that their origin in fact is same. 

These crossover phenomena within mean field approx- 



VII. APPENDIX: THE ELEMENTS OF THE 
MATRIX M 



In this appendix we give the expressions of the matrix 
elements (|29p in detail, 



M pp = p 2 x \C 2 M HH + 2CM HN (35) 
- 2C 2 M HN + (1-C) 2 M NN ), 

M CC = p 2 (M HH -2M HN + M NN ), (36) 
M pi = M ip =pCM Ht , (37) 
M Ci = M^ c = pM m , (38) 

M Cp = M pC = p 2 x (CM HH (39) 
+ (1 - 2C)M HN - (1 - C)M NN 

where the elements of the matrix are related to 



M HH = f3J mt p H ^ 



.d 2 k(x) 



dx 2 

1 a2 P 4t)-3i) 2 

l , ° P (l-r,)2 



2 Jmt£ 



dk(x) 
dx 



2 ' 



M 



NN 



Pph dp 2 H 

I . ° p (l-r,)2 



(3pN dpNdpN 



Mhn — Mnh 



dx 2 

Q2 p 4r] — 3r) 2 



p (l-r,)i 



dp H dp 



d 2 k{x) 
dxdx 



2JpH^^-+2J int p H Z, 
ox 



M N/i = M iN = 0. 



(41) 
(42) 

(43) 
(44) 

(45) 
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